





ae 








Institutional Archive of the Naval Postgraduate School 


Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1976 


Surface search radar tracking by a 
microcomputer Kalman filter. 


Wilson, Charles Howard 


http://hdl.handle.net/10945/17967 


Downloaded from NPS Archive: Calhoun 


| Calhoun is the Naval Postgraduate School's public access digital repository for 
D U DLEY research materials and institutional publications created by the NPS community. 
get Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS'‘s first 
KNOX appointed — and published — scholarly author. 


LIBRARY Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 






SURFACE SEARCH RADAR TRACKING BY A 
MICROCOMPUTER KALMAN FILTER 


Charles Howard Wilson 


Ly KNOX LIBRARY 
POSTGRADUATE SCHOVE 
EREY, CALIFORNIA 23948 





NAVAL POSTGRADUATE SCHOOL 


Monterey, California 





THESIS 


SURFACE SEARCH RADAR TRACKING BY A 
MICROCOMPUTER KALMAN FILTER 





by 


Charles Howard Wilson 


Lie) ewe 


Thesis Advisor: V. M. Powers 


Approved for public release; distribution unlimited. 





SECURITY CLASSIFICATION OF THIS PAGE (When Date Entered) 


| REPORT DOCUMENTATION PAGE 


| . REPORT NUMBER 2. GOVT ACCESSION NO 


4. TITLE (and Subtitie) 





READ INSTRUCTIONS 
BEFORE COMPLETING FORM 


3. RECIPIENT'S CATALOG NUMBER 






5S. TYPE OF REPORT & PERIOO COVERED 
Master’s Thesis: June 


1976 


PERFORMING ORG. REPORT NUMBER 















Surface Search Radar Tracking by a 
Microcomputer Kalman Filter 









8. CONTRACT OR GRANT NUMBER(e) 





7. AUTHOR(e) 









Charles Howard Wilson 






10. PROGRAM ELEMENT, PROJECT, TASK 
AREA & WORK UNIT NUMBERS 


12. REPORT DATE 
June 1976 


1S. SECURITY CLASS. (of thte report) 





: 9. PERFORMING CRGANIZATION NAME ANDO AODORESS 
Naval Postgraduate School 
Monterey, California 93940 









CON TROLLING OFFICE NAME ANO AOORESS 
Naval Postgraduate School 
Monterey, California 93940 






. MONITORING AGENCY NAME & AOORESS(If different from Controlling Office) 









Unclassified 


Se. OECLASSIFICATION 
SCHEOULE 


Naval Postgraduate School 
Monterey, California 93940 













COWNGRADING 


16. OISTRIBUTION STATEMENT (of thie Report) 


Approved for public release; distribution unlimited 


17. OISTRIBUTION STATEMENT (of the abstract antered in Block 20, if differant from Report) 


16. SUPPLEMENTARY NOTES 


19. KEY WORDS (Continue on reveree eide if necessary and identify by block number) 


Microprocessor, Microcomputer, Kalman Filter, Surface Search 
Radar, Radar Tracking 





20. ABSTRACT (Continue on reverse eide if neceseary and identity by bieck num 


U. S. Navy ships not equipped with NTDS currently perform sur- 
face radar target tracking manually, a tedious and inaccurate 


process. 
This thesis establishes that an 8-bit microprocessor, inter- 
faced to a common surface search radar, could compute the course 
and speed of multiple radar targets. The investigation includes 
the estimation of radar measurement errors, a determination of 
practical digital interface limitations, and the development 









DD ion, 1473 EDITION OF 1 NOV 88 1s OBSOLETE 
(Page 1) SECURITY CLASSIFICATION OF THIS PAGE (When Data Entered) 





a a cy 
Se cusMTy CLASSIFICATION OF THIS PAGESWhen Nete Entered. 





of a Kalman filter algorithm, a floating point arithmetic library 
for the INTEL 8080 microprocessor, and the microprocessor 
tracking system software. The error statistics and execution 
time of the microprocessor software are presented for several 
computer simulated target tracks. 


DD Form_ 1473 


7 an 4 
S/N 0102-014-6601 SECURITY CLASSIFICATION OF THIS PAGE(Mhen Data Entered) 





Surface Search Radar Tracking 


by a 
Microcomputer Kalman Filter 


by 


Charles Howard Wilson 
Lieutenant, United States Navy 
Enoeeunaversity of Idaho, 1970 


Submitted in partial fulfillment of the 
requirements for the degree of 


MASTER OF SCIENCE IN ELECTRICAL ENGINEERING 


from the 


NAVAL POSTGRADUATE SCHOOL 
June 1976 





NOX LIBRARY 
RADUATE SCHGuUG 
MUNTEREY, CALIFORNIA 23940 


ABSTRACT 


U. S. Navy ships not equipped with NTDS currently per- 
form surface radar target tracking manually, a tedious and 
inaccurate process. 

This thesis establishes that an 8-bit microprocessor, 
interfaced to a common surface search radar, could compute 
the course and speed of multiple radar targets. The 
investigation includes the estimation of radar measurement 
errors, a determination of practical digital interface 
limitations, and the development of a Kalman filter 
algorithm, a floating point arithmetic library for the 
INTEL 8080 microprocessor, and the microprocessor tracking 
system software. The error statistics and execution time 
of the microprocessor software are presented for several 


computer Simulated target tracks. 
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I. INTRODUCTION 


"In December 1973 INTEL Corporation shipped the first 
8-bit N-channel microprocessor, the 8080." /1_/ In the 
months that followed, several semiconductor manufacturers 
introduced their own version of INTEL's 8080 microprocessor. 
The resulting competition has driven microprocessor prices 
down and performance up. Advanced Micro Devices announced 
in April 1976, their AM9080A-4 microprocessor, with an 
instruction execution time one-half that of INTEL's 
original 8080, was available for $21.00 / 2_/. In the 
Same announcement, they predicted the AM9080A would cost 
$5.00 by 1980. Paralleling the advances in microprocessor 
technology, the manufacturers of Semiconductor memories are 
making steady progress in expanding the memory capability 
of single semiconductor chips, reducing memory access time 
and subsequently reducing the cost of high-speed semiconduc- 
tor memories used with microprocessors. 

These advances in semiconductor technology are taking 
place at a time when the U. S. Navy is debating the cost 
of rebuilding the fleet. A significant portion of the cost 
of building a modern warship is the electronics in its 
weapons systems. Microprocessors offer the potential to 
(1) reduce the cost of digital systems, (2) perform complex 
functions at remote stations, relieving the congestion at 
larger central computing facilities, and (3) perform func- 


tions currently handled by watch personnel, thus reducing 





the manning requirements of watch sections. 


Pee OSJBCTIVES OF THESIS 

The primary objective of this thesis is to demonstrate 
the capability of microprocessors to perform surface search 
radar tracking, thus relieving the tedium of this function 
from the typical underway watch section. 

The secondary objective of this thesis is to investigate 
the performance of a Kalman filter in the two-dimensional 
polar coordinate system, with measurement noise, typical 


of a surface search radar. 


Eee SCOPE OF THE WORK 

The scope of the work was limited to the development of 
a Kalman filter tracking algorithm for the INTEL 8080 
microprocessor and simulation of circuitry necessary to 
interface with the AN/SPS-10 surface search radar. 

A Kalman filter was chosen because it is computationally 
the most complex of the tracking filters typically used 
for radar tracking. Thus, it would be the most demanding 
on the microprocessor. Additionally, the Kalman filter is 
seldom implemented in its complete form due to its computa- 
tional complexity thus adding additional relevance to this 
investigation. 

The parameters for the AN/SPS-10 surface search radar 
were used for computations requiring specific radar para- 
meters. This radar was selected because it was available 


at the Naval Postgraduate School and is installed aboard 
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the majority of U. S. Naval ships. 

INTEL Corporation's 8080 microprocessor was used for 
the implementation of the Kalman filter because it was 
available at the Naval Postgraduate School along with sup- 
porting software. 

The interface between the microprocessor and the radar 
was not constructed due to: (1) Cost, time and manpower 
constraints, and (2) The performance of the system could 
not be effectively measured against non-cooperative targets 
me Opportunity. 

The performance of the radar and the interface equipment 
were effectively simulated on the Naval Postgraduate 
School's IBM-360 computer, thus permitting comparison of 
the microprocessor's performance against that of the IBM- 
360 system for known target tracks. 

The output of the tracking system would be displayed 
for human utilization rather than being electrically 
interfaced to fire control equipment. The minimum standards 
of performance were established as that obtainable from 
a well trained, fully manned underway watch section typical 
of that employed on a modern destroyer not equipped with 


a Naval Tactical Data System (NTDS). 


teal 





II. BACKGROUND 


A. PERFORMANCE OF A TYPICAL CIC WATCH TEAM 

Onboard naval ships not equipped with Naval Tactical 
Data Systems (NTDS), the tracking of radar contacts, both 
surface and airborne, is performed manually by the Combat 
Information Center (CIC) watch team. During normal peace- 
time steaming, this watch team may consist of from two to 
ten personnel, depending on the size of the ship and the 
complexity of its CIC equipment. Due to the risk of col- 
lision with other ships, their effort is primarily devoted 
to detecting, tracking and computing the course, speed, 
and closest point of approach (CPA) of other ships within 
a 20 to 30 nautical mile radius. This can be a tedious, 
error-prone task which generally requires a duplicate effort 
by the bridge watch team in order to avoid the awesome 
consequences of a collision at sea. The installation of 
an economical device to track and compute the tactical 
data for surface radar contacts could eliminate one or more 
watch-standers in CIC while improving the reliability of 
the information provided to the bridge watch team. 

The manual method of computing the course, speed and 
CPA of a radar target involves a graphical solution of the 
true and relative velocity vectors of the radar target's 
motion and own ship's course and speed. The target's 


relative position measurements are obtained from a radar 
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repeater. The range measurement accuracy is at best + 20 
yards with a well calibrated repeater and a well trained 
operator. The bearing measurement accuracy is at best 

+ 1 degree with a well trained operator. The own ship's 
course and speed used in determining the own ship's velocity 
vector are normally based on the ordered course and speed. 
The course actually steered by the helmsman may vary from 
the ordered course by + 2 degrees. Depending on the condi- 
tion of the ship's hull and its engineering plant, the 
speed will be accurate within + 1 knot. The combined 
results of these errors lead to a target course accurate 

to within + 10 degrees and a target speed accurate to 
within + 3 knots, not allowing for human error. fThe major 
disadvantages of the manual solution, in addition to its 
poor accuracy, is the time required to obtain the solution 
and the inability to obtain solutions if either the target 
or own ship maneuvers between radar position measurements, 
typically 3 or more minutes apart. The estimates of error 
are based on the author's personal observation of unalerted 


CIC watch sections. 


B. THE TRACKING PROBLEM 

The "Tracking Problem" as referred to in this investiga- 
tion is that of using radar measurements to accurately 
determine the course and speed of a waterborne radar target. 
The targets of interest include merchantmen, warships, and 


other ships capable of open ocean operations. As such, 
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their motion on the surface of the ocean can be approximated 
by a linear dynamic system. In a sampled data system, such 
aS a search radar where measurements are taken once each 
antenna rotation, the target dynamics can be represented 


in matrix form as / 3_/: 


Sn+1 : an id [ay 
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to n + 1 (E/W and N/S components) 


If the radar measurements were perfectly accurate and 


the target had zero acceleration, then the tracking problem 
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would be a simple velocity vector problem (Figure 1). 
Unfortunately, the radar bearing and range to the target on 
any particular scan is contaminated by errors which are a 
function of the radar system and the received signal to 
noise ratio. The measurement errors, to be discussed in 
the next section, are characteristically Gaussian. An 
optimal estimator is well suited to tracking problems such 
as this. 

"An optimal estimator is a computational algorithm that 
processes measurements to deduce a minimum error estimate 
of the state of a system by utilizing: knowledge of 
system and measurement dynamics, assumed statistics of 
system noises and measurement errors, and initial condition 
information.” / 4 / 

A Kalman filter is a common optimal filtering technique 
for estimating the state of a linear system /4/. It is 
chosen for implementation because it is computationally the 
most demanding as compared to other techniques used in 
Similar tracking systems and, as such, places the most 
demand upon the microprocessor. The details of the Kalman 
filter will be covered in Chapter IV following an analysis 
of the measurements and associated hardware which provide 


its inputs and determine its parameters. 


C. MEASUREMENTS 
In order to solve the tracking problem, it is necessary 
to make measurements of (1) target range, (2) target bearing, 


(3) own ship's course, (4) own ship's speed, and (5) time 
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of the measurement. It is assumed that all noise is Gaussian 
and received signal to noise ratios are large (>? 1). 

The target range is determined by measuring the time 
required for a radar pulse to travel to and from the 


target. 


R - Range to Target 
c - Speed of Light 
t - Elapsed Time from Pulse Transmission to 


Pulse Detection 
Using threshold detection to obtain a leading edge measure- 
ment, the theoretical standard deviation of error for the 
measurement of t is / 5_/. 

“i . A a 

2 ae : 
can?| = » W 
* (2 s/n)? 


Au = Error in Elapsed Time Measurement t 





——— 


~~. 


ig = Rise Time of Received Pulse 


s/n = Signal to Noise Ratio of Video Pulse 


If the rise time of the pulse Seales! is limited by the receiver 
IF bandwidth B, then ii ae = / 5_/. The theoretical 
standard deviation of range error using leading edge detec- 


tion is therefore: 
Cc 
— — 
°R ~ 2B(2 s/n)? 


OR = Standard Deviation of Range Error 
The mean range error is negligible when compared to the 


range measurement. Additionally, the mean error will tend 


inte 





to be cancelled by subtraction in the relative velocity 
calculations. Thus the computed target course and speed 
will not be affected by the mean range error. 

The target bearing is determined by measuring the bearing 
of the radar beam at the time the target video is detected. 
The mean error of the bearing measurement will be a function 
of the antenna beam pattern, and the target range and radar 
cross section. The resulting bearing to target will not 
be through the center of the target but the computed target 
course and speed will not be affected. 

Using beam splitting, a technique for estimating the 
center bearing of the target, the theoretical minimum stan- 
dard deviation of bearing error for a Gaussian beam distri- 
bution is / 5_/: 


1.06 @ 


2 2 =a 
ee Gc/n) =~ te 
2 ae: 
where 
Ty = standard Deviation of Bearing Error 
os = Two Way Beamwidth (= 0.4250 65, 8, = Half 


Power Beamwidth) (Same units as Tp) 
(s/n), = Signal to Noise Ratio at Center of Beam 
Ne = Number of pulses transmitted as the antenna 
rotates through 2 8, 
The standard deviation of both range error and bearing 
error are functions of the video signal to noise ratio. 
The video signal to noise ratio for a radar with single 


pulse detection is given by / 5_/: 
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eee -~°¢R/5 
P P x G rac 10 


n (baie at F_ BL, R 


woere 
Sie = Peak Transmitted Power (watts) 
G = Antenna Gain 
»\ = Wavelength (meters) 
o = Radar Cross Section of Target (m*) 
eo = Attenuation Constant of Propagation 


Medium (db/m) 
R = Range to Target (meters) 
kn ssoltzmanm-s Constant (1.23 x 107° joule/deg) 
T = Standard Temperature (290° K) 


O 
B,, = Receiver Noise Bandwidth (HZ) 
L, = System Losses 


When these relationships are applied to the AN/SPS-10 
(series) Surface Search Radar, with the AS-1161/SPS 
Antenna, the standard deviations for bearing and range can 
be related to target range and cross sectional area as 
Shown in Figures 2 and 3. The minimum performance charac- 


ters of this radar system are listed below. 


Peak Power = 190 kw L6/ 
Antenna Gain = 30 db Le 7 
Pulse Repetition Rate = 625 HZ (OTE 
Frequency = 5825 MHZ ow: 
Propagation Attenuation = .6 x 107? db/m LSS, 
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Noise Figure = 22.9 db Ley. 


Bandwidth = 5 MHZ La 
Noise Bandwidth (1.05 x B)= 5.25 MHZ L 5_/ 
System Losses = 5 db £5/ 
Antenna Rotation Rate = 17 RPM EE 
Antenna Horizontal Beamwidth = 1.9° Vea, 


As proposed earlier, a threshold detector could be used 
to indicate the occurrence of the leading edge of the 
received radar video pulse. The threshold to noise ratio 
of such a detector will establish the ranges at which 
various targets can be tracked and the minimum signal to 
noise ratio which the tracking system will encounter. 
Leading edge detection is preferred because it eliminates 
errors resulting from variations in the pulse width. 

When considering the threshold to noise ratio of such a 
detection scheme, it is necessary to establish the desired 
mean time between false alarms, and the desired minimum 
probability of detection of the target. 

The mean time between false alarms is a measure of 
how often the tracking system will erroneously recognize 
a noise pulse as a target. The tracking system can not 
function effectively if this occurs frequently. A mean 
time between false alarms of 60 sec has been selected. 

The mean time between false alarms is a function of the 


Signal to noise ratio and the IF bandwidth of the receiver 


5 _/. 
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ee 

y 2 

dk : ; 
(S45) = Threshold to Noise Power Ratio 

Bir = IF Bandwidth 

TEA = Mean Time Between False Alarms 


The exponential term in this expression is the inverse of 


the probability of false alarm (P,,) Vasey] 


(- Vn) eed 
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The mean time between false alarms as seen by the tracking 
system can be increased if the tracking system is disabled 


whenever targets are not anticipated. For the AN/SPS-10, 





the measurement circuitry could be disabled 77 percent of 
the time. This would multiply the mean time between 

false alarms by 4.32. An additional advantage can be 

Zained by uSing bearing and range gates to further limit 

the opportunity for false alarms. If such bearing and 

range gates isolated 1000 yard squares for each target 

and the tracking system were tracking eight targets, the 
mean time between false alarms would be increased by an 
additional factor of 1414. The combination of these factors 


leads to a probability of false alarm of 2.036 x 107? and 
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a threshold to noise power ratio of 10.8 or 10.3 db. 

The probability of detection of the target is a measure 
of how often the tracking system will recognize the target 
video when a target is present. It is essential to the 
effectiveness of the tracking system that it detect the 
target video reliably. The probability of detection is a 
complex function of the signal to noise ratio, the charac- 
teristics of the detector, and the probability of false 
alarm requiring numerical techniques / 5_/. Figure 2.7 
of Reference 5 shows a plot of the probability of detection 
of a sine wave in noise as a function of the signal to 
noise power ratio for various probabilities of false alarm. 
This plot shows that a signal to noise ratio of 12.5 db to 
14.0 db is required for a probability of detection of 0.90 
to 0.99, respectively, and a probability of false alarm 
of 2.0 x 107. If all targets were perfect radar reflec- 
fers, it would be possible to predict the performance of 
the threshold detector and the range at which various 
targets could be tracked. However, the reflected radar 
Signals from the targets of interest will fluctuate due to 
cancellation and reinforcement of the waves reflected 
from various portions of the target. These fluctuations 
can be modeled as either Swerling’s Case I or Swerling's 
Case II / 5_/. In either case, an additional signal to 
mease ratio of 8 db to 17 db will be required to insure 


a probability of detection of 0.9 to 0.99, respectively 


>_/ . 
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Thus the proposed threshold detection scheme should have 
a threshold to noise power ratio of 10.8 to achieve an 
effective mean time between false alarms of 60 sec. This 
would lead to a worse case minimum average signal to noise 
mero Of 20.5 db to 31.0 db for a probability of detection 
of 0.9 to 0.99, respectively. 

For a radar measurement system with the parameters 
established above, targets with the following average 
radar cross sections could be tracked with the AN/SPS-10 


at the ranges indicated. 


Average Radar Cross Range 
section (m@) (yards/meters ) 

il 3,600/3, 290 
i 6,400/5.850 
iLOP 11, 500/10, 520 
107 20,000/18, 290 
107 35 ,000/32,000 
10° 60,000/54, 860 


Based on a threshold to noise ratio of 10.8, the stan- 





(Figures 


2 and 3) wi respec- 





tively. 

The own ship's course and speed measurements come from 
the ship's gyro and pitometer log, respectively. Accordingly, 
their error statistics are a function of the equipment 
installed on the particular ship. This equipment is fre- 


quently associated with the ship's fire control equipment 
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and, as such, is intended for tracking operations. 

The gyro has a dual influence on the tracking problem. 
It not only determines own ship course but also determines 
the true bearing of the radar antenna. The mean error in 
radar bearing, and likewise own ship's course, will result 
in an error in the solution. This can be viewed as a 
shift in the reference plane which does not prevent the 
correct prediction of the future relative position of 
targets. As a refinement to the tracking system, it would 
be better to measure the average course and speed made 
good since the last measurement rather than the instan- 
taneous values at the time of the measurement. Such an 
average measurement could be achieved by continuous ly 
integrating the variations in the resolved N/S and E/W 
velocities using a time constant equal to one-half the 
period of the antenna rotation rate. This integration 
process would have the desirable side-effect of minimizing 


the variance of error in own ship's course and speed. 


For subsequent calculations, it is assumed that the standard 








The fifth and final measurement, the time of the measure- 
ment, could be made with a very high accuracy, depending on 
the quality of the clock used. The objective is to deter- 
mine the elapsed time between position measurements. The 
antenna rotation rate determines the mean elapsed time 


between measurements on a unique target, and the pulse 
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repetition rate of the radar determines the useful precision 
to which the elapsed time needs to be measured. The mean 
rotation rate for the AS-1161/SPS Antenna is 17 RPM. The 
mean pulse repetition rate for the AN/SPS-10 (series) 

Radar is 637.5 HZ. Thus, the mean elapsed time between 
measurements is 3.53 sec and the useful precision of the 


time measurement is + 0.784 m sec. 


D. MICROPROCESSORS 

A microprocessor is a large scale integrated (LSI) 
semiconductor device or family of devices which contain 
an arithmetic logic unit (ALU) and several working registers 
capable of performing binary or BCD arithmetic and/or 
Boolean logic operations. Microprocessors generally fall 
into one of five categories: (1) single-chip controllers, 
Meet-bpit family, (3) 8-bit family, (4) 16-bit family, 
and (5) bit slices / 9_/. 

The single chip controllers are LSI devices which 
incorporate all the necessary hardware, including program 
and data memory circuitry, on a single chip. Their data 
processing applications include devices such as calculators, 
gas pumps, cash registers, and scales. They include devices 
with 4-bit and 8-bit ALUs. Typically the 4-bit versions 
perform serial BCD operations. They are limited in their 
Speed of program execution by the power of their instruc- 
tion sets and the size of their ALUs. The complexity of 


the operations they can perform is limited by the size of 
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their memories. Some manufacturers are developing single- 
chip versions of their 8-bit family devices which will 
expand the capabilities of the single-chip controllers 
meethe point of overlapping the 8-bit families / 9_/. 

The 4-bit microprocessor families are similar in 
function to the single-chip controllers except the memory is 
not included on the chip. Therefore, the 4-bit families 
are capable of handling larger programs. The term family 
refers to the set of compatible chips which go together 
to make up the microprocessor. They require greater design 
effort than the single chip controllers because several 
chips, including clocks and I/O interface chips, along 
with the separate memory, are required for even the simplest 
task. Like the 4-bit controllers, the 4-bit microprocessor 
families are limited in the speed of their program execution 
by the power of their instruction sets and the size of 
their ALUs. 

The 8-bit families are presently dominating the micro- 
processor applications / 9 /. Their instruction sets are 
approaching the complexity of those used in minicomputers. 
They typically provide for 8-bit data manipulation along 
with 16-bit memory addressing. "---INTEL Corporation's 
8080 family, with its enhanced 8080A CPU, Motorola Semi- 
conductor's 6800 family, with its enhanced 6800D CPU, 
and Rockwell's PPS-8 family currently rank one, two, and 
three in popularity among users." / 9 7 


The strength of the 8080 lies in the power of its 
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instruction set and the extent of its software support. 

The 8080 has 72 basic instructions, including conditional 
branching, decimal as well as binary arithmetic, logical, 
register to register, stack control and memory reference 
instructions /1_/. These instructions can be executed 

in 4 to 18 cycles of its 2 MHZ clock allowing some instruc- 
tions to be executed in 2 usec /1_/. These instructions 
facilitate 16-bit addition and register to register trans- 
fers. Among its software support are a System Monitor, an 
Assembler, a Text Editor, and a PL/M compiler, cross 
assembler, and simulator. The popularity of PL/M, an 
INTEL derivative of PL/1, and the 8080 have caused other 
manufacturers to model their microprocessors after the 
8080, such as Advanced Micro-Devices 9080A, or to develop 
their own version of PL/M such as PL/M 6800 by Intermetrics 
for the Motorola's 6800. 

The 8-bit devices are limited in performance by their 
8-bit data path to and from memory. This necessitates 
using several memory cycles to read or write address and 
other 16-bit variables. 

The 16-bit microprocessor family, including National's 
PACE, Texas Instruments Inc.'s 9900, and General Instrument 
Corp.'s CP1600, attempt to reduce the number of instructions 
and machine cycles required to manipulate memory addresses 
(typically 16-bit variables), and data originating from 
A-D converters or going to D-A converters (typically 10-16 


bits). Their performance approaches that of minicomputers. 
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Western Corp.'s MCP 1600 was originally designed for DEC's 
LSI-11 16-bit minicomputer and was subsequently distributed 
as an independent device / 9_/. 

The bit-slice architectures, Such as Advanced Micro 
Devices Inc.'s AM 2900, INTEL Corp.'s 3002, and Texas 
Instrument Inc.'s SBPO400, are taking on the tough process- 
control and high speed controller jobs previously handled 
by minicomputers. A typical bit-slice device is a 2 or 4 
bit vertical slice of a CPU which is horizontally expandable 
meee or 4 bit increments to form a 16 bit, 32 bit or other 
mcr ple of 2 or 4+ bit CPU. They utilize Schottky TTL, 


tf 


L or ECL circuitry and thus are much faster than the 
previously mentioned microprocessors which are primarily 
MOS devices. The bit-slice architecture offers the ulti- 
mate in design flexibility while demanding the most effort 
by the designer. The design effort may include an extensive 
amount of microprogramming prior to development of the . 
working program. Additionally, bit-slice devices are the 
most expensive form of the microprocessor. 

In summary, the 8-bit and 16-bit microprocessors are 
the best alternatives to the tracking problem. The single- 
chip controllers and 4-bit families have insufficient memory 


and/or are too slow to handle the Kalman filter algorithm. 


The bit slice architectures have excess capability while 


requiring unnecessary design and hardware cost. As will 
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PTE.) SYorEM INTERFACES 


The required interface equipment can be grouped into 
three categories; (1) the radar interface, (2) the own 
ship course/speed interface, and (3) the operator/user 
interfaces. Although the interface hardware was not con- 
structed, it is essential to the validity of the proposed 
tracking system to verify the feasibility of making the 
required radar and own ship measurements within the pro- 


posed tolerances. 


A. RADAR INTERFACE 

The radar interface serves four functions; (1) to con- 
vert the radar video to rectangular pulses compatible 
with the logic circuitry in the range measurement subsystem, 
(2) to identify the radar video of the target whose position 
is to be measured, (3) to measure and store the bearing, 
range and time of measurement of the designated radar 
target, and (4) to transmit a signal to the own ship 


course/speed and the operator/user interfaces to indicate 


when the radar measurements have occurred. 
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voltage would be a technician adjustment made by first 


setting the IF amplifier gain to obtain a good video 
display on the radar's master repeater and then adjusting 


the threshold voltage on the comparator to just detect a 
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aud ae accomplished either 
with software or hardware. software solution would result 
in the measurement and processing of a large number of 
targets which the operator had no interest in and would 


place an unnecessary burden on the microprocessor. Thus 





The width and time of occurrence of the bearing and 


range gates can be determined in several ways, all of which 
are a function of the accuracy with which the future loca- 
tion of the particular target can be estimated. This 
estimate is a function of (1) the maneuverability of the 
target, (2) the elapsed time between measurements and (3) 
the accuracy of the computed course and speed of the target. 
The majority of targets of interest will be capable of 
turning radii of 100-1000 yards, velocities of 0-40 knots, 


and accelerations of 0-.5 yards/sec*. The antenna rotation 
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rate determines the mean time between measurements, which 


in the case of the AS-1161/SPS antenna equates to 3.53 


seconds £ 8_/. The accuracy of the computed course and 








Because of the uncertainty of these many factors, it is 
necessary to use wide bearing and range gates to insure 
that the target is not lost. The unfortunate outcome of 
using wide bearing and range gates is the increased likeli- 
hood that two adjacent targets will enter each others 
bearing and range gates and become confused. 

The best solution to this dilemma is to allow the 
operator to select a target classification from his control 
station which would determine which of several predeter- 
mined bearing and range gate widths or width determining 


algorithms would be used with the designated target. isom 
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digitizing process is carried out with the digital zero 
equal to the analog zero it will introduce a mean error 
equal to one-half of the value represented by the least- 
significant bit and a standard deviation of error equal to 
372 times the value represented by the least significant 
bit. If the digital zero is offset by minus the previous 
mean error, then the mean error will be zero and the 
standard deviation of error equal to 1272 times the value 


represented by the least significant bit. Obviously, the 


zero offset approach is preferred. 





For the system under consideration, the least signifi- 


cant bit should represent the standard deviation of error 





and the useful precision of the time measurement (+ 0.784 


msec), as a minimum, The maximum value of each measurement — 


for the time measurement. If an 8-bit microprocessor is 


to be used, then data words of 8, 16, or 24 bit lengths 





are practical. In addition, the INTEL 8080 microprocessor 


has some instructions for manipulating 16-bit data words. 
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Thus 16-bit data words were considered optimum. This 


results in the following data format: 


Measurement Range MSB LSB 
ering O-27f rad tY rad 0.0959 mrad 
Range 0-60,000 yds 30,000 yds 0.9155 yds 
Time 0-8.0 sec 4.0 sec 0.1221 msec 


By changing the most significant bit to the next even power 


of two, the units can be rationalized as follows: 


Measurement Range MSB LSB 
Bearing 0-277 rad 4.0 rad g7t3 rad 
Range 0-65,535 yds Bem oo vas 1.0 yds 
Time 0-7.999 sec 4.0 sec 2713 sec 


The resolution resulting from this format will result in 


less than 1.0 percent increase in the standard deviation 


Smeerror as a result of digitizing the data. The remaining 





surement could be conveniently made 
by combining one of the synchro-to-digital converter modules 
such as Transmagnetics series 1651, which has 14-bit 
resolution, and two 16-bit latchs, enabled by the video 
Circuitry, to store the first and last bearings from 
which video was detected from a particular target on 
each scan. The microprocessor software could then compute 


the center bearing to the target. An alternate scheme 


oY, 





would be to use a Single latch to store the first bearing 
from which video was detected from a particular target on 
each scan and a counter to count the number of video pulses 
from the target on each scan. The microprocessor software 
could then use this data to estimate the center bearing. 

A third and far less accurate scheme would be to simply 
use the bearing of the first video pulse. The first 


scheme is recommended. 





range measurement is essentially a time interval 
measurement. A range of 1 yard corresponds to 6.1 nsec. 
Thus a binary counter would have to be clocked at 1L(6)3} ane, 
MHZ to count radar range in intervals of 1 yard. Ordinary 
Transistor-Transistor-Logic (TTL) counters are limited to 
clock frequencies of 30 MHZ and below, while high-speed TTL 
counters can operate as high as 50 MHZ /10/. Schottky- 
Clamped TTL counters can operate reliably at clock fre- 
quencies of 100 MHZ my eeGteten—Coupled-hogzic (ECL) 
Seecultey could perform at higher clock frequencies but 
the complexity and cost of the circuitry would increase 
considerably due to the unavailability of ECL Medium-Scale 
Integrated (MSI) semiconductor chips with the desired 
functions. A Shottky-Clamped TTL Counter, such as the 
74S197, is suggested, with a compromise reduction in the 
clock frequency to 81.96 MHZ resulting in a minimum range 
resolution of 2 yards. An alternate scheme would be to 
mix ECL circuitry with the Shottky counters. 


Another consideration in the range measurement counter 


38 





is the propagation of carry pulses through the counter. The 
74S197 binary counter has a maximum propagation time of 13 
nsec through each stage. For the 15 stages required for a 
Maximum range of 65,536 yards, a total propagation time of 
195 nsec would be required for a carry propagation to the 
highest significant bit. This is 16 times the clock period, 
which would lead to unstable data at the counter outputs. 
The easiest solution to this dilemma is to stop the clock 
input to the counter at the time the target video is 
detected and store the counter value after it has stabilized. 
This would necessitate a separate range counter in the 


interface hardware for each target being tracked. 





easurement Miele: is comparatively straight 
forward. A 16-bit counter driven by a 8,192 HZ clock, with 
a storage latch for each target will suffice. As in the 
range counter, the carry propagation delay must be taken 
into account when storing the time measurement. Due to 
the slower clock frequency, the counter output will be 
stable 99 percent of the time. Simple combinational logic 
circuitry could be employed to synchronize the time measure- 
ment with the clock to avoid storing unstable counter 
outputs. 

Using currently available hardware, the radar interface 
can be accomplished with resolutions well within the 


radar measurement errors specified earlier in this section. 


The following data formats are recommended. 


ay 





Measurement Range hicnmot Bats Resolution 


Bearing 0-27 rad 14 e395 mrad 
Range 0-65,535 yards 15 2.0 yards 
Time 0-8 sec 16 21221 msec 


B. OWN SHIP COURSE AND SPEED INTERFACE 

The own ship course and speed interface serves two 
functions; (1) to smooth the variations in own ship's 
course and speed, and (2) to measure and store the smoothed 
course and speed to be associated with a particular target. 
There are two options of how to perform these functions, 
stemming from the fact that the microprocessor software 
uses own ship's course and speed data solely to compute 
own ship north/south and east/west velocities. 

The first approach would be to process the course and 
speed measurements independently. The course measurement 
could be smoothed either mechanically or with an over-damped 
servo feedback system. A synchro to digital angle converter 
could be employed to digitize the measurement. The speed 
measurement could be smoothed electrically and digitized 
uSing an analog to digital converter. 

The second approach would be to combine the course and 
speed measurements into north/south and east/west velo- 
cities. The course synchro signal could be converted to 
two analog signals corresponding to the sine and cosine of 
own ship's course using a synchro to DC sine/cosine converter, 


such as Transmagnetics model 655N V (series). These analog 


Lo 





Signals could then be multiplied by an analog speed signal 
with two DC analog multipliers. lLow-pass filters could be 
used to smooth the two signals and then digital to analog 
converters used to digitize them. This approach would 
Simplify the work of the microprocessor but additional 
errors would be generated in the analog multipliers. 

A compromise approach would be to use the first system 
but substitue a synchro to digital sine/cosine converter 
in place of the digital angle converter. With this tech- 
nique, the microprocessor would receive smoothed north/ 
south and east/west course component measurements and 
Smoothed speed measurements. 

Assuming a standard deviation of error in the course 
measurement of 0.1 degree (1.745 x 107? radians) prior 
to digitizing, a digital data format of 11.8 bits would be 
required for minimum resolution. Using a data format 
similar to that used for the target bearing (MSB = 2.0 
radians) and 14-bit resolution, the standard deviation 
of error will only increase 0.3 percent as a result of 
digitizing. Assuming a standard deviation of error in the 
speed measurement of 0.1 knot prior to digitizing, and a 
maximum speed of 100 knots, a digital data format of 10 
bits would be required. Using a 12-bit analog to digital 
converter or a 12-bit synchro to digital converter with the 
most significant bit representing 64 knots, the standard 
deviation of error will increase only 1.6 percent, as a 


Result of digitizing. 
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In summary, the own ship course and speed interface 
could be accomplished using available hardware, with the 


following data formats and measurement precision. 


Measurement Range No. of Bits Resolution 
Course 0-27 rad 14 OG lpmlbad 
Speed 0-128 knots LZ .03125 knots 


C. OPERATOR/USER INTERFACE 

The operator/user interface serves two functions; (1) 
to provide operator control of the tracking system, and 
2eco display the output data to the users. 

The operator should have control over (1) which radar 
video is to be tracked, (2) the number of targets being 
tracked, and (3) the tracking parameters to be used with 
each target. 

It is proposed that the operator would control the 
system in the following manner. The measurement interface 
Circuitry would be divided into tracking elements, each 
capable of handling one target. Each element would be 
composed of a bearing and range gate generation circuit, 

a range measurement circuit, and storage circuits for the 
bearing, range, time, own ship's course and own ship's 
Speed data. 

The operator's control panel would have an off/on/start 
Switch for each tracking element. The system would be 
turned on with all tracking elements turned off. To track 


a target, the operator would place his radar repeater's 
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bearing and range cursor on the target video to be tracked 
and momentarily place the tracking element switch in the 
start position when the radar sweep passed the target and 
then place the tracking element switch in the on position. 
This process would momentarily couple the radar repeater's 
bearing and range cursor video to the tracking elements 
bearing and range gate circuit and enable the tracker's 
measurement circuitry, thus taking a measurement of the 
target's position. If the measurement circuitry detected 
carmeeet Video, an indicator light would be lit. To stop 
tracking a target, the operator would simply turn the 
tracking element off. A control data word (8 bits) would 
be periodically read by the microprocessor to advise it as 
to which tracking elements have updated measurement data 
available. The operator would have a second multi-position 
Switch which would allow him to select the bearing and 
range gate widths and the variance of acceleration para- 
meters to be used with each target. 

The output data to be displayed would include (1) target 
designation, (2) target bearing, (3) target range, (4) tar- 
get course, and (5) target speed, as a minimum. For more. 
sophisticated systems the target's closest point of approach 
(CPA) could be computed and displayed as (1) CPA bearing, 
(2) CPA range and (3) time of CPA. The speed of the data 
output is critical to the capability of the system to track 
multiple targets. Direct memory access techniques are 


suggested as the only practical solution. 
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After the display data is extracted from the tracking 
system, several display techniques could be used. Tabular 
displays in the form of 7-segment LED modules or cathode 
ray tubes would facilitate displaying the data to several 
users at the operators station and also transmission of 
the data to remote stations, such as the bridge, CO's 
cabin, and weapons stations. One of the major advantages 
of an automatic tracking system is its capability to 
instantly and continuously disseminate tactical data to 
remote stations thus eliminating the need for passing 
this information on sound powered phones as is currently 


the practice on non-NTDS ships. 
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IV. KALMAN FILTER 


As stated earlier, the Kalman filter is a common optimal 
filtering technique for estimating the state of a linear 
system. "Among the presumed advantages of this type of 


data processor are that it minimizes the estimation error 











in a well defined statistical sense and that it utilizes 


— a 





all measurement data plus prior knowledge about the system. 
The corresponding potential disadvantages are its sensitiv- 
ity to erroneous a priori models and statistics, and the 


inherent computational burden." Ll 


A. DISCRETE KALMAN FILTER ALGORITHM 

The rotation of the surface search radar antenna 
results in periodic measurements of the target positions. 
Accordingly, the discrete form of the Kalman filter, rather 
than the continuous form, is appropriate. The discrete 
Kalman filter can be described by a set of matrix 
equations. 

Using the target state cae) representation presented in 
the description of the tracking problem, in Section I B, 


the matrix form of the measurement process is / 3_/: 


w. = HZ +v 
Nn Nn Nn 
x, (n) 
W == 
eal) 
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ae Vale = Measured x Coordinate at Scann 








y,, (1) = Measured y Coordinate at Scan n 
aie 0 -407,0 
2 aia 
Z., = Target State Vector at Scann 
vn “festa 
ee ve a 
ve 
v,(n) = Random Measurement Error on x(n) 


vy (n) = Random Measurement Error on y, (n) 


The estimated state vector at scan n (Z_) As ae 


7 =6? 
an _ Zn-1 


$ - (See Tracking Problem Page 14) 


A 
Z 


i a Optimal Estimate of Target State Vector 


at Scan n-l 


The optimal estimate of the target state at scann (Z.) 


Heesiven by / 3_/: 


Z =2 + K( HZ) 
an Kal n'“n 7 Zn 
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The gain matrix (K_) is computed from / 3_/: 


4 ES coll poe al 
K = Pt (HP_H 


+ R y7i 
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Estimated Covariance of Estimation Errors 


ack 
T 


nN Matrix Prior to Processing Scan n Measurements 
2 2 
o,. (n) Ree (n) 
R = 3 5 = Covariance of Measurement 
Soe (n) oy (n) Error Matrix at Scan n 


Following the processing of scan n measurements, an 
optimal estimate of the covariance of estimation error 


matrix (PB) for scan n can be computed / 3_/: 


“A 
Ia 


nee ee K HP 
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Identity Matrix 


From this, it is possible to compute the estimate of 


miemecovariance of estimation error matrix for scann?-tti 


ae . 


= Poel ao wT Cen 
earl agree yy errr) n+l Taare 


Q = 0° = Variance of Random Target Acceleration 


Mies basis for the discrete Kalman filter in this form 


is the legitimacy of the target maneuver model. The 





is true for the targets of interest. Thus, this form of 
the discrete Kalman filter is valid if the position 
measurements were made in the x and y coordinates. 

Using the polar coordinate system shown in Figure 4, 
the following transformations can be used with the radar's 


bearing (6) and range measurements (R). 
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North 


® (Bearing) 
R (Range) 


West 


SoOuGch 


FIGURE 4 COORDINATE SYSTEM 
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x R sin @ 


m 


IH 


ve R cos @ 


Using these transformations, the elements of the covariance 


of measurement error matrix are / 11_/: 


2 2 2 


Zoe 2a 
oe = oF Simmcma: i %% cos 8 
ee Zz 2 Z DE 
oy = or BOS ae le % Ssin™ @ 
2, Z 2 Z Cee 
ie ise - fr Sis $+ sin 2 0 
07." = Variance of Range Error 
a = Variance of Bearing Error 


This modification completes the description of the discrete 
Kalman filter. 

In some applications, it 1s possible to precompute the 
values of the gain matrix and store them in the computers 
memory for subsequent use. This greatly simplifies the 
software algorithm. There are two reasons why this is not 
possible in this case. 

The first is that the elapsed time between measurements 
(T) is not constant. This is primarily due to the fact that 
the radar antenna rotates relative to own ship's heading. 
Thus, changes in own ship's course can significantly 
increase or decrease the elapsed time between measurements 
2) 


The second reason is that the covariances of measurement 
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error (OL, Fei ae are a function of the bearing and range 
to the target. Consequently, the gain matrix must be com- 
puted on line for each target being tracked. A flow chart 
of these calculations is shown in Figure 5. 

The solution of the Kalman filter equations in matrix 
form is unnecessarily complex for the microprocessor. 
Therefore, an algebraic solution for the value of each 


element in the various matrices was formulated. Appendix 


A lists this algebraic form of the Kalman filter equations. 


B. KALMAN FILTER PARAMETERS 

The performance of the Kalman filter is determined by 
three parameters; (1) variance of bearing error (7°), 
(2) variance of range error orn and (3) variance of 
target acceleration acy. The first two parameters 
are functions of the radar, the ship's gyro, and the inter- 
face system. The third parameter is a function of the 
target's relative motion and the desired smoothness of the 
output data. 

The variance of bearing error is equal to the sum of 
the variance of the radar bearing error, (.0058 deg”), the 


variance of the bearing- interface digitizing error (.00016 


dee“), and the variance of the ship's gyro error (0.01 deg“). 





The variance of range error is equal to the sum of the 


Z 


variance of the radar range error (50 yd™ or 42 m°) and 
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the variance of the range interface digitizing error 


er. 3 ya’ or lat m). 


error is equal to 51 ya© (43m). | 


The variance of target acceleration parameter is a 





compromise between the user's desired smoothness of the 
output data for non-maneuvering targets and the ability of 
the filter to provide the closest estimate of a target's 
position and velocity for a maneuvering target. 


A possible solution to this dilemma is to allow the 





the bearing and range gate widths, could be selected by 


enabling one of several read-only-memories (ROMs) containing 
the desired parameters. The particular ROMs installed 

at any one time could be selected from a library of such 
ROMs by the Commanding Officer depending on the ship's 
mission. 

The determination of what values of variance of accelera- 
tion to be used with various targets are best determined by 
modeling LUA A computer model of the tracking problem 
was written and used to investigate the system performance 
with various values of variance of acceleration. The 
results of this investigation are presented in the figures 


discussed in Section IV E. 


C. KALMAN FILTER INITIAL CONDITIONS 


As illustrated by Figure 5, there are two basic loops in 
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the Kalman filter algorithm. They represent the recursive 
calculations of the optimal estimates of the covariance of 
Bermation error matrix (B) and the target state matrix (Z). 
An initial condition must be established for each of these 
based on a priori information of the measurement process and 
the targets of interest. 

The initial condition of the target state matrix is 
equal to the expected value of the initial state. The 
expected values of the initial position elements, Xo and 


Yo» are the first position measurements: 


Xo 


R(O) sin 9(0) 


aig R(0) cos @(0) 


@) 
The expected values of the initial velocity elements, Xo 
and Yo» are zero. Thus, the filtering of the target state 
does not begin until the second measurement. 

The initial condition of the covariance of estimation 


error matrix is equal to the covariance of error of the 


initial state estimate. 


Aw ae -_ ww == 7 
Pe 202, - 24/7 L225 - 26S 
BY, )-_/ 0 0 0 
mn 0 EL (X,)* 7 0 0 
= een 
0 0 0 Ey Meare) 7/ 0 
0 0 0 ES (%,)94 


The non-zero elements in this matrix are functions of 


the initial relative velocities which might be encountered, 
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and the measurement errors of the initial position. Assuming 
a uniform probability distribution of relative velocities 
from -50 yd/sec to +50 yd/sec (+ 46 m/sec), and an initial 
range of 50,000 yds (46,000 m), the variance of initial 
Peeition error is 12,000 yas* CHO OO0 m“) and the variance 
of initial velocity error is 830 ya“/sec* (700 m“/sec“) 
‘ for both the X and Y coordinates. 

The primary influence of the initial conditions is 
the time required for the filter to achieve a steady-state 
solution. Large values lead to longer settling times for 
nonmaneuvering targets where small values lead to longer 


settling times for maneuvering targets. 


D. SYSTEM COMPUTER MODEL 

In order to evaluate the performance of the Kalman filter 
with various targets and parameters, a computer model of 
the system was derived. The computer model included the 
Simulation of target and own ship motion, the radar measure- 
ment and interface process, the Kalman filter, and the 
computation of the filtered target course and speed. 
Additionally, the program included computation of measure- 


ment errors, target track errors, and associated statistics. 





motion by digitally integrating the north/south and east/ 
west velocity components using an integration step of 0.01 
sec and double precision arithmetic. Target and own ship 


course changes were accomplished by a second-order nonlinear 
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control system model with the ordered course as the input, 
a limit on the maximum angular rate, and the actual course 
steered as the output. The control system model parameters 
and angular rate limit were adjusted to simulate typical 
course changes for the targets of interest. Speed changes 
were accomplished by a first-order linear control system 
model with the ordered speed as the input and the actual 
speed as the output. No attempt was made to duplicate 
actual ship maneuvers, but only to approximate the maneuvers 


typical of the targets of interest. 





3; were modeled by first 
comparing the rotation of the radar antenna to the bear- 
ing to each of the targets. Upon detection of a target, 
the bearing and range to the target were noted along with 
own ship's course and speed, and the time of detection. 


Random numbers representing the Gaussian distribution of 


measurement errors were added to the bearing and range 


measurement. The measurements were then truncated to 


Simulate the digitizing process of the interface system. 





implemented through FORTRAN 
statements Similar to the algebraic expressions listed in 
Appendix A. The only modifications were those related to 


improving the efficiency of executing the computations. 





from the relative north/south and east/west positions 
and velocities provided by the Kalman filter, and the 


simulated measurement of own ship's course and speed. 
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E. KALMAN FILTER PERFORMANCE 


Using the computer model previously described, the per- 
formance of the Kalman filter was evaluated for two classes 
of targets, maneuvering and nonmaneuvering. In all cases, 
own ship was proceeding on a steady course (000°) and at 
a fixed speed of 20 knots (37 km/h) from an initial position 
at the origin of the cartesian coordinate system. Eight 
targets were initially positioned and proceeding on the 
courses and at the speeds shown in Table 1. A unique set 
of random numbers was used with each target for the various 
computer runs. This eliminated the influence of variations 
in the random number distributions on the statistics of 
the Kalman filter output. 


The following system parameters were used: 


RADAR PARAMETERS 


Pulse Repetition Rate OO. 0 pps 
Range Precision 7.0574 yds 
(6.453 m) 
Bearing Precision 0.1258 deg 
Antenna Rotation Rate WO (Ol Tegan 


DIGITAL INTERFACE PRECISION PARAMETERS 


Range LSB 4.0 yds 
(3.6576 m) 
Bearing LSB 0.02197 deg 
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Time LSB 0.1221 msec 
Own Ship Course LSB 0.02197 deg 


Own Ship Speed LSB 0.03125 knot 
(0.05788 km/h) 


The Kalman filter measurement error parameters (ca and 


o-.*) were selected to match the computed measurement errors 


resulting from the radar and interface parameters. 


one =ib.o70 x 107° rad’ (0.01599 deg“) 
o.* = 51.14 ya* (42.76 m®) 


The Kalman filter initial conditions were the same as those 
discussed in Section C of this chapter. The variance of 
target acceleration (07%) parameter was chosen as the inde- 
pendent variable. 

Both scenarios began with: the acquisition process for 

each of the eight targets, which involved setting the 
initial conditions of the Kalman Hilcer. ands ran tor. 5 
minutes of simulated time. For the nonmaneuvering case, 
no changes were made throughout the run. For the maneuver- 
ing case, each of the targets executed a 45 deg course 
change beginning 120 sec after the start of the run. 
Targets 1, 4, 6 and 8 turned right, while targets 2, 3, 
5 and 7 turned left. These turns were executed with an 
angular rate of 3.0 deg/sec. Figure 6 shows an example 
of the plot developed during such a left turn. 

The actual measurement error statistics, based on 680 


measurement samples from the computer model, for the 
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nonmaneuvering and the maneuvering target scenarios are 
listed in Table II. There was good correlation between the 
measured variances of error and the computed variances. 
Additionally, the mean errors were approximately equal to 
minus one-half the values of the least significant bit 
used in the interface Simulation. This is in agreement 
with the theoretical mean error generated by truncation. 
Each scenario was repeated for seventeen logarithmically 
distributed values of the Kalman filter variance of target 
acceleration parameter from 0.0005 ie fee (0.00042 m“/ 
Bec *) ico ot nO fae (0.85 WC /oe In all cases, the 
Kalman filter achieved a steady state condition in 60 sec, 
corresponding to 17 target position measurements. Steady 
state is defined as the condition where the Kalman filter 
gain matrix (K,) would be constant for targets with constant 
relative position and constant sample intervals. The 
bearing, range, course and speed of each target was com- 
puted from the Kalman filter's optimal estimate of target 
state and compared with the actual target bearing, range, 
course and speed. Error averages were computed for the 
time interval from 60 sec to 300 sec. Figures 7, 8, 9 
and 10 show plots of the average root mean Square of error 
for the target bearings, ranges, courses and speeds, 
respectively, as a function of the variance of target 
acceleration parameter. It should be noted that the 
error averages for the maneuvering target are a function of 


the elapsed time between maneuvers and the time spent in 
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maneuver. Thus, the plots for the maneuvering target only 
apply for the particular scenario used. 

Several conclusions can be drawn from these figures. 
For nonmaneuvering targets, the variance of acceleration 
parameter (07) should be small (< .001 yd-/sec* (.0008 
fyeec’)) if smooth course and speed values are desired. 
For maneuvering targets, there is an optiaml value of 
o* e210 sl / ee” (.08 e/sech below which the filter 
does not follow the target maneuvers well and above which 
the target data will be noisy regardless of the target 
maneuvers. There is a lower limit to how well the Kalman 
filter will track a maneuvering target even with an optimal 
value of i Finally, the bearing and range to a non- 
maneuvering target will be improved by the Kalman filter, 
as shown by the standard deviation of measurement error 


line on Figures 7 and 8, over a wide range of on, * values. 
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V. MICROPROCESSOR PROGRAM 


A MCS-8 microcomputer system, using INTEL’s 8080B 
microprocessor, was used to implement the Kalman filter 
algorithm. The objective of implementing the algorithm 
was to determine limits to the time required to accomplish 
the Kalman filter computations and the accuracy of the 


results. 


A. PROGRAMMING LANGUAGE 

INTEL Corporation has developed the PL/M compiler, a 
derivative of PL/1, for use with their 8080 microprocessor. 
PL/M is a block structured language which speeds micro- 
processor program development by relieving the programmer 
of the tedious task of memory and register management (bl Je 
It facilitates the use of variable names while still permit- 
ting the programmer the freedom to reference specific 
memory locations. The basic instruction set includes 
arithmetic and Boolean algebra operation on 8-bit and 16-bit 
variables and/or constants. Program control is accomplished 
with Procedure, Do, Go To, If, Call and Return statements. 
Logic functions such as comparison, shift and rotation are 
also included. Floating point and transcendental functions 
must be implemented by subroutines. PL/M permits the 
program control provided by assembly languages while 


taking a much more readable form. 
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B. FLOATING POINT ARITHMETIC 

Floating point arithmetic is required whenever the range 
of numbers represented by the program variables is large or 
unpredictable / 12_/. Such is the case with the Kalman 
filter algorithm. This is particularly true if the operator 
is given the freedom to alter the filter parameters. A 
disadvantage of using floating point arithmetic is the 
additional execution time required. 

A suitable set of floating point arithmetic procedures 
were not available. The procedures listed in Appendix B 
were derived from several existing procedures and the 
general guidelines detailed in Reference 12. 

A 16-bit mantissa, 7-bit exponent and 1-bit sign were 
chosen to represent the floating point numbers using three 


8-bit data words in the following format: 


‘Mantissa (16-bits) . Sign (1-bit) 
Exponent (7-bits) 


M - Memory Address Associated with the Floating Point 
Number 
The 16-bit mantissa was chosen because it is the largest 
binary number that INTEL's 8080 microprocessor can manipu- 
late without a considerable increase in the complexity of 
the progr4mming. Additionally, it is compatible with the 
input and output data. The mantissa is always left justi- 


fied with the binary point located to the left of the most 
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Significant bit. Thus the mantissa is always equal to or 
greater than 0.5, and less than 1.0. 

The sign was located in the third word in order that the 
mantissa be a full 16-bits. A positive sign is represented 
by a binary 0 and a negative sign is represented by a 
binary 1. 

The remaining 7-bits of the third word are used for the 
exponent. The value of the exponent is equal to 64 plus the 
power of two to be associated with the mantissa. 

The data and constants in the PL/M programs are written 
in this format with hexadecimal numbers. The following 


are examples of this form. 


DECIMAL HEXADECIMAL FLOATING POINT 
Ons 80H, OOH, 4OH 
jee 80H, OOH, 41H 
| ae 80H, OOH, C1H 
10.0 AOH, OOH, 44H 
64.0 80H, OOH, 47H 
O.1 CCH, CDH, 3DH 


The largest magnitude number that can be represented 
is FFH, FFH, 7FH, which corresponds to 9.223 x 1018, The 
smallest magnitude number that can be represented is 80H, 


20) Zero is 


OOH, OOH which corresponds to 2.710 x 10° 
represented by OOH, OOH, OOH. 
The typical execution time of the arithmetic and trig- 


onometric procedures listed in Appendix B, using INTEL's 
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8080B microprocessor, are listed below. 


FLOATING POINT PROCEDURE MAXIMUM EXECUTION TIME (msec) 


Addition (ADD) 5S 
euotcraction (SUB) 1.4 
Multiplication (MULT) 3.4 
Division (DIV) 4.1 
Comparison (COMPARE) OniS 
Square Root (SQRT) nL 50. 
Cos (TRIG 1) 18.0 
Sin (TRIG 2) HS. © 
Cos and Sin (TRIG 3) 300 € 
Arctangent (TRIG 4) 19.0 


The actual execution times of these procedures are a func- 
tion of the numbers involved. 

Overflow and underflow provisions are included in the 
arithmetic procedures. Overflow will result in a solution 
equal to the largest possible number with the appropriate 
Sign. Underflow will result in a solution equal to zero. 

The teen of the sine and cosine functions, using 
the procedure "TRIG", is at worst + 0.000488 if the argu- 
ment is in the range 0 to 277. The precision of the arc- 
tangent function, using the procedure "TRIG", is equal 
to + 0.0610 mrad. 

The four arithmetic procedures are used by passing the 
variables’ addresses to the procedure with a call statement. 


The following PL/M statements are examples. 
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FUNCTION PL/M STATEMENTS 


met Y = 2 Came ADE er, «0% 2Z): 
on = Z OeulL SUN (any aes, AA 
hex B= C Coli ApoB, «C): 
ii = 0 Can tei (a. <M, <Q): 


Note: The dot symbol proceeding a variable name in PL/M 
indicates the address of the variable. Some 
variables may correspond to addresses of other 


Variables, in which case the dot would be omitted. 


The "COMPARE" procedure returns a binary integer quantity 


which indicates the relationship of the arguments. 


PL/M STATEMENT 
N = COMPARE (.A, .B); 


RESULT 

If A<B Then N = 0 
If A=bB Then N = 1 
If A>B Then N = 2 


The argument for "SQRT" must be positive. 


FUNC TION PL/M STATEMENT 
y=Jfx Call SQRT (.X, .Y); 


The procedure "TRIG" can be used to compute the sine 
and/or cosine, and arctangent by passing the addresses of 
the arguments plus a function indicator flag to the pro- 


cedure with a call statement. It requires less time to 
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compute Sine and cosine of the same angle with one call 
rather than two. The following are examples of the multiple 


uses of the "TRIG" procedure. 


FUNCTION PL/M STATEMENTS 

fe COSmGA) Coit nen Ge. W054, 2.) 
YS Sim A CotemerrG MCaW ey «Ay 2) 3 
m= Cos (A) and 

= Sala (04) Cami eTRECM(ok =e YaeeeA,. 3) 
A = Arctan (Y/X) Calle TREGHG@ xs oA, 4b); 


Note: W is a dummy variable used to fill the unused argu- 
ment location in the call statement when only cosine 


Se sine isedesired. 


C. TRACKING SYSTEM PROGRAM 

An experimental version of the tracking system software 
was written to determine the execution time of the Kalman 
filter algorithm and the error statistics of the resulting 
target track data. This program, written in PL/M, is 
listed in Appendix C. There are three functional parts 
femithis program; (1) data input, (2) the Kalman filter, and 
(3) computation of output data. 

The data input portion of the program includes the 
Simulated interface system output data, and the simulated 
control panel functions. The simulated interface system 
data was generated in the correct floating point format 
by the FORTRAN program used to model the system performance. 
The variables "BRG," "RNG," and "TIM" correspond to the 
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measured target bearing, range and time of detection, 
respectively. The own ship north/south and east/west 
velocity are represented by the constants "NSVEL" and 


(a 2 Z 


"EWVEL." The Kalman filter parameters Oe» and C,. 


ao 
are represented by "SIGA," "SIGTH," and "SIGR." The 
variables "MARKOLD" and "CONTROL" serve the functions of 
distinguishing which target data corresponds to existing 
tracks, and which control panel switches are turned on. 
These variables are associated with the variables "SIMFLAG 
( )" and "DFLAG" which correspond to the interface unit 
from which the target data came and thus the target's track 
number. In an actual operating system the program would 
use INPUT statements or direct memory access to read in the 
target measurement and control panel switch data. The 
measurement data would have to be converted from fixed 
point to floating point format before processing. 

The Kalman filter statements are grouped into two long 
DO statements. The first is used to initialize the filter 
when a new target is first processed. The second group is 
used to process new data for existing target tracks. The 
sequence of statements and variable names parallel the 
algebraic expressions listed in Appendix A. Temporary 
variables were introduced where necessary to avoid repeti- 
tive operations and to store intermediate results. 

Using the Kalman filter optimal estimates of target 
position and relative velocity, the third portion of the 


program computes and displays the target bearing, range, 
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course and speed. 

Following the program debugging, a series of tests, 
using a full range of typical input data, were run to 
verify the correct computation of each variable. The 
FORTRAN model was used to compute the "correct" values for 
each variable using the Naval Postgraduate School IBM 360 
computer. It was found that the computational errors were 
confined to a range of significance sufficiently small to 
avoid instabilities in the filter. After processing several 
measurements, the Kalman filter estimates of target state 
Peputed with the microprocessor were within the same 
range of error as the IBM 360 solutions. The covariance 
of estimation error and Kalman gain matrix elements computed 
by the microprocessor converged rapidly on the values 
computed by the IBM 360 computer. 

As a check on the microprocessor's performance a series 
of eight test runs were made, each with a different target. 
The eight simulated targets were identical to the non- 
maneuvering targets used previously except the signs of 
the east/west initial positions for targets 3-8 were 
negative. This change was made because the execution time 
and precision of the "TRIG" procedure is worst for targets 
in the north-west quadrant. Nonmaneuvering targets were 
used because the available memory space limited the amount 
of data that could be loaded with the program. The test 
runs were restricted to 130 sec tracks. A variance of 


acceleration of .00O1 aC ae ( .00084 i feces was selected 
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to insure that computational errors would not be over- 
shadowed by normal estimation error in the output. It 
was assumed that these tests would represent worst case 
eenadi_tions. 

The test program occupied a total of 2E00, ¢ bytes of 


memory. This was divided as follows: 


Floating Point Procedures 98A, ¢ Bytes 
Kalman Filter Algorithm ig ae Bytes 
Bearing, Range, Course and 

Speed Computations 2204¢ Bytes 
Variable Storage 36346 Bytes 


The error averages for the test data are listed in 
Table IIIl. These error averages represent the performance 
of the microprocessor tracking program after the Kalman 
filter achieved a steady state condition, approximately 
50 sec after initialization. This resulted in a relatively 
small sample set of 20-21 samples for each target. The 
error statistics of the input data during this time are 
listed in Table IV. As compared to the measurement error 
statistics presented in Table II, there is less correlation 
with the theoretical measurement errors for the microprocessor 
test data. Therefore, the error averages listed in Table 
III should be considered typical only for a short period 
of operation rather than for long periods with all possible 
targets. 

An analysis of these results reveals the microprocessor's 


influence on the track data. By comparing the course and 
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speed data for targets 3, 5 and 7 with the data for targets 
4, 6 and 8 it can be seen that for targets with similar 
course and speed the errors increase with range. Addition- 
ally, errors are larger for targets with small relative 
velocities. The combination of long range and small rela- 
tive velocities can be expected to result in larger errors 
due to the limitations of the floating point arithmetic to 
handle large numbers with small differences. 

The final parameter of interest is the execution time 
of the microprocessor program. The INTEL 8080B micro- 
processor was able to execute 85 consecutive iterations of 
Gieekalman filter algorithm in 37.5 secs, resulting Tay etch al 
average execution time of 0.44 secs. The computation of 
target course and speed required a worst case execution 
time of 54 msec. It would not be necessary to compute 
the target course and speed for every target on each antenna 
scan. The additional processing required to output the 
course and speed data would best be accomplished by a 
separate microprocessor dedicated to each display device. 
In general, a network of microprocessors could be used to 
share the task and thus increase the number of targets 


which could be tracked. 
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VI. CONCLUSION 


The tedious, error-prone task of tracking surface search 
radar targets manually, as is common onboard naval ships 
not equipped with NTDS, can be eliminated with the use of 
large scale integrated microprocessors. 

The AN/SPS-10 Surface Search Radar, installed on most 
U.S. Navy ships, could provide target position measurements 
Watn Sufficient accuracy to perform automatic tracking of 
surface radar targets. Using leading edge detection and 
currently available hardware it would be possible to 
digitally measure the bearing and range to a target with 
a standard deviation of error of 0.126 deg and 7.15 yards 
(6.54 m) respectively. 

Using a 4-state Kalman filter and appropriate conver- 
Sions for polar coordinate measurements, it 1S possible to 
obtain optimal estimates of target north/south and east/ 
west position and velocity. From this data and knowledge 
of own ship's motion, target course and speed can be com- 
puted far more accurately than by manual tracking tech- 
niques. 

A set of floating point arithmetic procedures was 
developed for INTEL's 8080 microprocessor based on a 
floating point number format with a 16-bit mantissa and 
a 7-bit exponent. These procedures maintain sufficient 


accuracy to reliably compute the Kalman filter variables 
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while limiting the execution time sufficiently to permit 
complex functions, like the Kalman filter, to be implemented. 

The Kalman filter, and Reus and speed calculations were 
successfully implemented using the microprocessor develop- 
ment system at the Naval Postgraduate School. - The average 
execution time of the Kalman filter was 0.44 secs while 
the course and speed calculations required 54 msec. This 
was accomplished with computed course and speed errors for 
nonmaneuvering targets typically less than 2.0 deg and 0.5 
knots. With the 3.5 sec sweep of the AN/SPS-10, a single 
system could continuously provide course and speed data 
for seven targets. This could be increased by using faster 
microprocessors currently entering the market. 

It is recommended that future work in this area be 
devoted to (1) improvements in the floating-point software, 
2) development of the interface hardware, and (3) improve- 
ments to the Kalman filter algorithm. The floating point 
software could be improved by reducing the execution time 
of the arithmetic procedures and increasing the precision 
of the transcendental functions. The Kalman filter algorithm 
could be improved by minimizing the effects of own ship's 
motion, and reducing the execution time. 

A microprocessor-based tracking system has the potential 
to provide any ship with a digital automatic surface track- 
ing capability, without the expense of NTDS size equipment. 
If implemented, it could reduce the manning of underway 


watches while improving the quality of the tactical 
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information available to the Officer of the Deck. There 
are additional applications of similar systems for such 


functions as sonar tracking and navigation. 
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